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ABSTRACT 

We investigate the deceleration of Solar Energetic Particles (SEPs) during their propagation 
from the Sun through interplanetary space, in the presence of weak to strong scattering in a 
Parker spiral configuration, using relativistic full orbit test particle simulations. The calculations 
retain all three spatial variables describing particles’ trajectories, allowing to model any transport 
across the magnetic field. Large energy change is shown to occur for protons, due to the combined 
effect of standard adiabatic deceleration and a significant contribution from particle drift in the 
direction opposite to that of the solar wind electric field. The latter drift-induced deceleration is 
found to have a stronger effect for SEP energies than for galactic cosmic rays. The kinetic energy 
of protons injected at 1 MeV is found to be reduced by between 35 and 90% after four days, and 
for protons injected at 100 MeV by between 20 and 55%. The overall degree of deceleration is 
a weak function of the scattering mean free path, showing that, although adiabatic deceleration 
plays a role, a large contribution is due to particle drift. Current SEP transport models are found 
to account for drift-induced deceleration in an approximate way and their accuracy will need to 
be assessed in future work. 

Subject headings: Solar Energetic Particles, drift, deceleration 


1. Introduction 

Solar Energetic Particles (SEPs) are com¬ 
monly measured near 1 AU in the interplanetary 
medium, in most cases at large distances from the 
region where they were accelerated, in solar flares 
and Coronal Mass Ejection (CME) driven shocks. 
Understanding their propagation through space is 
important to gain information about the accelera¬ 
tion processes and to test theories of charged parti¬ 
cle transport in turbulent magnetic fields, relevant 
to a variety of astrophysical environments. Since 
SEPs are a radiation hazard for satellite infras¬ 
tructure and humans in space, accurate modelling 
of their propagation is important also to develop 
a Space Weather forecasting capability. 

During their motion through interplanetary 
space, SEPs undergo changes in their kinetic en¬ 
ergy. At the present time, a well established 
framework for the description of kinetic energy 
change of SEPs exists, based on the theory of adia¬ 
batic deceleration within a transport equation, ini¬ 


tially developed to describe galactic cosmic rays. 
In its classic description, adiabatic deceleration 
results from the cooling of an isotropic energetic 
particle gas as it moves o ut wards in the ex panding 
solar wind ( Parkeilfl965lb ISkillind ( 197lh derived 
a focussed transport equation including adiabatic 
deceleration for the case when the plasma con¬ 
vection velocity is at an angle with respect to 
the m agnetic field, retaining pitch-angle as a vari¬ 
able. iRuffolol (1995) obtained an alternative for¬ 
mulation for a Parker spiral geometry, again con¬ 
sidering particles of arbitrary pitch-angle. Both 
approaches use a deceleration term within a trans¬ 
port equati on. Th e term s derived by [Skilling 


(Il97llf and iRuffolol (Il995il are equivalent for a 
Parker spiral magnetic field. 

The majority of modern models of SEP prop¬ 
agation include _ ad iaba tic deceleration _by 


us- 


i ng e ither the IRuffolol ([1995) or the ISkillina 


(|l97ll) formali sm within t h e foc us sed transport 


equation [e.g. lLario et al.l (119981) : lAgueda et al 
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(12008 ): Zhang et al. (12009 ): Droee et al. ( 2010l) : 
Mason et alJ ( 20121) : Wang et al.1 ( 20121) ]. Numer¬ 


ical solution of this equation shows that adiabatic 
deceleration needs to be taken into account when 
studying SEP propagation, as it affects the decay 
phase of events sign ifican tly f or pr oton energies up 
to about 20 MeV ( Ruffoloifl99a ). Intensity pro¬ 
files are affected in a complex way tha t varies in 


the d iff erent stages o f the e vent (IKocharov et al 


1998). Mason et al. ( 20121) modelled the trans¬ 


port of heavy ions and showed that the effects of 
adiabatic cooling are very strong below 1 MeV 
nucleon -1 . They pointed out that, as a result 
of deceleration, particles measured within a given 
energy channel of a detector located at 1 AU may 
have had a much higher energy when originally 
injected at the Sun. 

The standard focussed transport approaches 
( Skilling! 1971 : Ruffololll995 ) assume that energetic 
particles remain tied to magnetic field lines dur¬ 
ing their propagation, so that only a single spatial 
variable, the distance travelled along a magnetic 
field line, needs to be retained to describe the par¬ 
ticle’s trajectory. This is a common assumption in 
many studies involving SEP modelling and data 
analysis. 

However in two recent publications we pointed 
out that drifts due to the gradient and curvature 
of the Par ker spiral magnetic field are importan t 
for SEPs ( Marsh et al. I2 oI 3 iDalla et al.1 2013 ). 


This result is based on full orbit test particle sim- 
ulations of SEPs in a Parker spiral configuration 
( Marsh et, al. ( 20131 ). Paper I) as well a s on analyt¬ 
i cal ex pressions for the particle drifts ( Dalla et al.l 
( 2013 ). Paper II). As a result of drift, particles ex¬ 
perience transport in longitude and latitude, per¬ 
pendicular to the magnetic field line on which they 
were originally injected. The scattering condi¬ 
tions, described by a mean free path A, do not 
strongly affect the amount of drift experienced. 
In fact, significant curvature drift is found to be 
present even when particles propagate scatter-free. 


In this paper we examine an important con¬ 
sequence of drifts for SEPs: the strong particle 
deceleration due to drift motion in the direction 
opposite to that of the solar wind electric field, 
studied here for the case of protons. The electric 
field is present in the fixed (observer) frame due to 
the outward propagation of the solar wind plasma. 
Drift-induced deceleration has been studied previ¬ 


ously in t he context o f galactic co smic rays (Kota 
19791 : Barnden fc Bercovitchl Il97 5;). but has not 


been discussed for SEPs up to the present time. 
We show that, unlike adiabatic deceleration, drift- 
induced deceleration is present also in the absence 
of scattering. We evaluate whether the stan dard 
focus s ed transport equation formalisms ([Skilling 
1971 : R.uffolol 1995) describe the effects of drift- 


induced deceleration on the particle distribution 
function. 


2. Energy change for SEPs as derived from 
test particle simulations 

In Paper I, the propagation of SEPs in a Parker 
spiral configuration in the presence of scattering 
was studied by means of a relativistic full orbit 
test particle code. The output of the simula¬ 
tions showed that, rather than remaining confined 
within the magnetic field lines delimiting the in¬ 
jection region at the Sun, particles experience sig¬ 
nificant drift (see e.g. Figures 3 and 5 of Paper I). 
The amount of drift, measured in a local Parker 
spiral coordinate system as the displacement from 
the magnetic field line through the initial posi¬ 
tion, increases with energy, i.e. it is much larger 
for 100 MeV than for 1 MeV protons. It is par¬ 
ticularly large for partially ionised heavy ions due 
to proportionality of grad B and curvature drift 
velocities to the charge to mass ratio. 

Here we examine the results of simulations sim¬ 
ilar to those presented in Paper I, focussing on the 
change in kinetic energy experienced by the par¬ 
ticles during their propagation. We follow four 
initially monoenergetic populations, each consist¬ 
ing of 1000 protons, characterised by values of the 
initial kinetic energy Kq of 1 MeV, 10 MeV, 100 
MeV and 1 GeV. The particles are injected from 
an 8° x 8° region at 1 solar radius and centered at 
latitude 5=20°. Their initial velocity directions 
are randomly chosen in a semi-hemisphere in ve¬ 
locity space pointing away from the Sun. Particles 
are scattered with an average frequency defined by 
a scattering mean free path A. During each scat¬ 
tering event the particle’s velocity is randomly re¬ 
assigned from an isotropic distribution. Particles 
are scattered in the solar wind frame and in this 
frame energy is assumed to be conserved during 
a scattering event. This is a reasonable approxi¬ 
mation since the majority of scattering events in 
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our simulations take place far from the Sun. For 
comparison, scatter-free runs are also performed. 
Other parameters of the simulations have been de¬ 
scribed in Paper I. 

The Parker spiral magnetic field used in the 
simulations is unipolar (pointing outwards from 
the Sun, as was assumed in the simulations of Pa¬ 
per I), i.e. the presence of the heliospheric cur¬ 
rent sheet is not modelled. The Lorentz equation 
of motion for each particle is solved in the fixed 
frame, in which the Sun is rotating with angular 
velocity fi and the solar wind is moving radially 
outward at constant speed v sw . This is the frame 
which is relevant for comparison of simulation re¬ 
sults with observations taken by spacecraft. 

Figure [l] shows the kinetic energy change expe¬ 
rienced by particles of the four initially monoen- 
ergetic populations (one for each panel, with ini¬ 
tial energy Kq as labelled) at time f/=102 hours, 
plotted as a function of the change in colatitude 
experienced by each particle. The mean free path 
is A=1 AU at all energies. Kinetic energy change 
AK=K-K 0 is shown as percentage of the initial 
energy Kq , as a function of A9=9—9q where 9 is 
the particle’s colatitude and 9q its initial value. 
Particles in all four populations are seen to have 
experienced strong deceleration over the time of 
the simulation, apart from a minority of those with 
I\q =1 GeV, whose kinetic energy increases. 

It should be noted that the scale on the pr¬ 
axes of the panels of Figure Q] is different for the 
four populations, showing that the 1 MeV parti¬ 
cles have a much smaller A 9 compared to 1 GeV 
particles. The observed change in colatitude is 
caused by drift due to the curvature and gradient 
of the Parker spiral magnetic field (Paper I, Paper 
II), and the amount of drift increases with particle 
energy. 

Figure [T| shows a trend for the deceleration to 
increase with A 9 for the cases where Kq= 1 MeV, 
10 MeV and 100 MeV, but not for the case Kq= 
1 GeV. 

In Figure [2] the kinetic energy change averaged 
over each of the four particle populations, is plot¬ 
ted versus time. The top panel shows the case A=1 
AU, corresponding to the same simulations shown 
in Figure [Q while the bottom panel presents the 
time dependence of energy change in the scatter- 
free case. Figure [2] shows that the decrease in ki- 
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Fig. 1.— Percentage change in kinetic energy as a 
function of colatitude change A 9, four days after 
injection, for protons of four initially monoener- 
getic populations with initial energy Kq indicated 
in each panel. Particles are injected from an 8° x 8° 
region at latitude 5= 20° and the mean free path 
is A=1 AU. 
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Fig. 2.- Time variation of AK/Kq averaged over 
each initially monoenergetic population with Kq= 
1 MeV, 10 MeV, 100 MeV and 1 GeV. The top 
panel is for A=1 AU and the bottom panel for the 
scatter-free case. Other parameters of the simula¬ 
tions are as in Figure [l] 


netic energy is fast during the first ^20 hours of 
the propagation, and continues to take place more 
slowly at later times. Deceleration is seen also in 
the absence of scattering, although it is of a re¬ 
duced magnitude compared to the A=1 AU scat¬ 
tering case. 


3. Deceleration processes 


3.1. Energy change due to drift in the elec¬ 
tric field potential 


As energetic particles propagate through the in¬ 
terplanetary medium, in the fixed, non-rotating 
reference frame they are subject to the electric 
field associated with the outwardly moving solar 
wind. The solar wind electric field is given by: 


E = —- 


x B 


(1) 


where v sw is the solar wind velocity, B is the inter¬ 
planetary magnetic field (IMF) and c is the speed 
of light. 


We consider B to be a unipolar Parker spi¬ 
ral, and the solar wind flow to be radial, uniform 
over the Sun’s source surface and time indepen¬ 
dent. With these assumptions, the e xpress ion for 
the electric field takes t he form (Burns & Halpern 


1968; Winge & Coleman 1968): 


A r 2 

E =-O Bq — sin 9 eg (2) 

c r 

where eg is the standard unit vector in a (r, 9 , <j >) 
spherical coordinate system centered on the Sun, 
with r the radial distance, 9 the colatitude and 
(j> the longitude. Here Bq is the magnitude of the 
magnetic field at a reference radial distance r$ and 
S2=2.87xl0~ 6 rad s -1 is the solar rotation rate. 
A= 1 when B points outward from the Sun and 
and A=—l when B points inwards. 

The potential associated with the electric field 
of Eq. © is given by: 

A 

$ =-HBofn cos 9 (3) 

c 


In the absence of other effects causing deceler¬ 
ation, the kinetic energy change due to motion in 
the direction of the solar wind electric field is given 
by: 

AA"($) = DB 0 rQ (cos 9 — cos0 o ) (4) 
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with q the particle charge, 6q its initial colati¬ 
tude and the subscript (<f>) indicating that this 
is the energy change associated with the potential 
of Eq. © . Expressions equivalent to Eq. © were 
derived bv iKotal ( 1979h . 


If particles do not move in colatitude, 0=9q and 
the kinetic energy change given by Eq. 0 is zero. 
However Papers I and II demonstrated that curva¬ 
ture and grad -B drift combine to produce a signif¬ 
icant component of drift in 6 1 which, from Eq. 0, 
will cause kinetic energy change. The drift in co¬ 
latitude is present also in the absence of scattering. 

Eq. 0 can be rewritten as: 


AK ($) = A sign(g) £ 0 (cos 6 — cos 6 0 ) (5) 

where: 

£o = MfiB 0 rg ( 6 ) 

For a proton or an electron, using the same values 
for Bq and ro used in Papers I and II, we obtain 
that £o—247 MeV. 

Paper II showed that for the case of a positive 
polarity of the Parker spiral IMF (A>0) and par¬ 
ticle initial positions not far from the heliographic 
equatorial plane, ions drift downwards (increas¬ 
ing colatitude) and electrons upwards (decreasing 
colatitude). Hence from Eq. @ both ions and 
electrons will experience a negative i.e. a 

deceleration, as a result of drift motion. 

For a negative polarity configuration (A<0) 
protons near the equatorial plane drift upwards 
and electrons downwards, so that the effect of the 
drift is again deceleration for ions and electrons. 

The expression for the energy change (Eq. ©) 
does not involve the initial particle energy explic¬ 
itly, however the latter is an important parameter 
in determining how much the particle drifts in co¬ 
latitude, i.e. the value of cosd — cos do- 


3.2. Energy change due adiabatic deceler¬ 
ation 

In standard transport theory, heliospheric en¬ 
ergetic particles experience cooling due to the ex¬ 
pansion of the solar wind, which carries the tur¬ 
bulence that produces particle scattering. Formal¬ 
ism to describe this adiabatic deceleration within a 


transport equation was derived by Parker (1965) 
for the case of an isotropic particle distribution 
function and radially directed magnetic field. 


Particle scattering in the solar wind frame is 
essential for adiabatic deceleration. In the absence 
of scattering, this process does not take place. 


Ruffolol ( 1995h generalised the classic treatment 


to the case of particles of arbitrary pitch angle and 
a Parker spiral configuration, pointing out how 
the combination of focussing and scattering pro¬ 
duces a deceleration in the solar wind frame. The 
rate o f deceleration for a single particle is given by 
( Ruffolo 1995 1: 


(V 

dt 


= ~P v s 


sec ip 


2 L(z 


(1 — u' 2 ) + cos ip— sec ibu‘ 
dr 


(7) 

where pf and p! are the particle’s momentum and 
pitch angle in the solar wind frame, ip is the an¬ 
gle between the Parker spiral magnetic field line 
and the radial direction and L(z) is the focussing 
length. Here z is the only spatial variable retained 
within the model, describing the distance travelled 
along a magnetic field line. 

From Eq. 0, the characteristic time of decel¬ 
eration T——{l/p'dp’/dt)~ l depends only on the 
geometrical parameters, including the colatitude, 
via the angle ip and L, and on the particle pitch 
angle, i.e. it does not depend explicitly on the par¬ 
ticle’s kinetic energy nor species. However the lat¬ 
ter parameters influence the overall trajectory and 
amount of time that a particle spends at a given 
location, and as a result the overall energy change 
experienced. If the scattering mean free path is 
small, particles tend to spend a longer time in re¬ 
gions near the Sun, where adiabatic deceleration 
is stronger, so that a lower value of A gives rise 
to more dec eleration associated with this process 
( Parker 19651 ). 

An adiabatic deceleration formalism alt ernati ve 
to Eq . 0 had been derived earlier by ISkillind 
(119711) . for a generic magnetic field. For the case of 
a Parke r spi ral fi eld, it can be shown analytically 
that the Skil 
form as the 

formulations are equivalent in this case, as noted 
in the latter paper. 


igj (El97 


) exp ression takes the same 
Ruffolo (jl995l ) one, so that the two 


4. Interpretation of simulation results 

Having introduced the concepts of drift-induced 
and adiabatic deceleration we can now interpret 
the results shown in Figures [T] and [5] Our simu¬ 
lations do include adiabatic deceleration because 


5 










































A=0.3 AU A= 1 AU A=10 AU 


0 



0 

\ 

0 


1 MeV ' 

o 

CM 

1 

(%) 



o 

CM 

1 

(%) 

\ 

o 

CM 

1 

(%) 


' 

^ -40 



s 2 - 40 

4i‘ 

-40 







a v 




^ -60 


1 

^ -60 

afl r 

c 

co 

1 


- 

< 



< 


<1 



-80 



-80 

* 

-80 


- 










-1.0-0.5 0.0 0.5 1.0 1.5 2.0 -1.0-0.5 0.0 0.5 1.0 1.5 2.0 -1.0-0.5 0.0 0.5 1.0 1.5 2.0 

Ad (deg) Ad (deg) Ad (deg) 



-101254 -101254 -101254 

Ad (deg) Ad (deg) Ad (deg) 



0 5 10 15 20 25 0 5 10 15 20 25 0 5 10 15 20 25 

Ad (deg) Ad (deg) Ad (deg) 



0 20 40 60 80 100 120 0 20 40 60 80 100 120 0 20 40 60 80 100 120 

Ad (deg) Ad (deg) Ad (deg) 


Fig. 3.- Percentage change in kinetic energy as a function of Ad , four days after injection, for three values 
of the scattering mean free path A and varying initial kinetic energy. The red dots are from simulations 
including both adiabatic deceleration and drift-induced deceleration, the green dots are for simulations 
without adiabatic deceleration, i.e. scattering taking place in the fixed frame. The blue square represents 
the scatter-free case. 
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scattering is taking place in the solar wind frame. 
They also include drifts which produce changes in 
colatitude and thus energy change. 

In Figure [0 the observed trend for decelera¬ 
tion to increase with Ad (as seen for Kq= 1 MeV, 
10 MeV and 100 MeV) is consistent with drift- 
induced deceleration playing a significant role. 
The fact that deceleration is present also in the 
scatter-free case (Figure [5] bottom panel), when 
no adiabatic deceleration is taking place, confirms 
this statement. 

Low energy SEPs experience very little drift, 
however their initial kinetic energy is much smaller 
than £ 0 (given by Eq. ©), the magnitude of the 
largest possible deceleration associated with drift 
processes. Considering Eq.©, one can see that 
even a small A 8 can produce a significant percent¬ 
age change in their energy, the behaviour shown in 
Figure ID (e.g. for Kq=1 MeV). SEPs of initial en¬ 
ergy 1 GeV typically drift by tens of degrees, how¬ 
ever their drift-induced deceleration as percentage 
of initial energy is much less than for lower en¬ 
ergy particles due to the larger ratio Kq /£o ■ As 
one considers particles of even higher energy, drift- 
induced deceleration becomes less significant. 

The deceleration seen in Figures Q] and [2] is a 
combination of adiabatic and drift-induced decel¬ 
eration. However it is possible to ‘switch off’ adi¬ 
abatic deceleration by scattering particles in the 
fixed frame, rather than in the solar wind frame. 
When this is done, particle energy is conserved in 
the fixed frame during a scattering event, while a 
change in pitch angle does result, influencing the 
particle’s trajectory. Scatter-free simulations rep¬ 
resent an extreme case: here no adiabatic deceler¬ 
ation takes place due to the absence of scattering 
events and particles propagate together as a beam. 
Drift and the deceleration it produces cannot be 
eliminated in our simulations. 

Figure [3] shows the output of simulations car¬ 
ried out in three very different scattering condi¬ 
tions, each corresponding to a column in the fig¬ 
ure. Red points are the results for the cases with 
adiabatic deceleration (scattering carried out in 
the solar wind frame). Points in the central col¬ 
umn (A=l AU) are the same as in Figure [1] and 
the columns on the left and right show the re¬ 
sults for A=0.3 AU and 10 AU respectively. The 
light green points show the energy change ob¬ 
served when scattering is carried out in the fixed 


frame, i.e. no adiabatic deceleration is occuring, 
for the same values of the mean free path. The 
blue squares represents the deceleration observed 
in the scatter-free case: here all particles experi¬ 
ence the same A 8 due to drift, resulting in a single 
point in the plots and no adiabatic deceleration. 

Considering first of all the red points in Figure 
0 one can see that as the level of scattering in¬ 
creases (i.e. moving from right to left in the figure), 
more deceleration is observed and the range of per¬ 
centage change in AA' increases. The amount of 
drift, represented by the range of A 8 values also 
increases, as was noted in Paper I. Considering 
that in our simulations the mean free path values 
vary by almost 2 orders of magnitude, from 10 AU 
to 0.3 AU, the amount of deceleration appears to 
be only a weak function of the mean free path. 

When no adiabatic deceleration is present (light 
green points), the amount of deceleration observed 
has much less scatter and follows a line which is 
found to coincide exactly with that predicted by 
Eq. ©, i.e. the observed energy change is purely 
drift-induced deceleration. It is also apparent that 
any effects due to large finite Larmor radius are 
small in our simulations, as the data points for 
this case lie along a single line in Figure [3] without 
significant scatter. 

We can now compare the results for the cases 
with adiabatic deceleration, without it and the 
scatter-free situation. Looking at the low scatter¬ 
ing regime (A=10 AU), it is apparent that the ad¬ 
dition of adiabatic deceleration broadens the dis¬ 
tribution around the line predicted by Eq. ©, and 
the effect becomes more marked for stronger scat¬ 
tering. Switching on adiabatic deceleration elimi¬ 
nates the low values of |AAT/Ao| which are present 
in the case of purely drift-induced deceleration 
(e.g. for AT 0 =1 MeV values of |AA7A' o |<30°). 
In the strong scattering case, adiabatic deceler¬ 
ation is more efficient, resulting in a distribution 
of AK/Kq values which does not show a clear de¬ 
pendence on Ad, unlike for the A=10 AU case. 
The weak dependence of the overall deceleration 
on the value of the mean free path demonstrates 
that drift-induced deceleration plays a very signif¬ 
icant role. 

Regarding the reason why a small minority of 
particles with A'o=l GeV show acceleration (Fig¬ 
ure [U bottom panel), we note that it is possible 
for the process responsible for adiabatic deceler- 
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ation to produce acceleration in the fixed frame 
in some cases (as can be seen from the qualita¬ 
tive d escription accompanying Figure 1 of Ruff'olo 

(H995ll V 


In the simulation discussed so far, particles were 
injected at 1 solar radius. We have also carried out 
simulations with injection at 40 solar radii for A=1 
AU and the deceleration observed is very similar 
to that shown in Figure [T] Therefore the effects 
describe above apply both to particles energised 
in solar flares and to those accelerated in CME 
driven-shocks at distances up to at least 40 solar 
radii. 


5. Discussion and conclusions 

Our relativistic full orbit test particle simula¬ 
tions of proton propagation in the Parker spiral 
have allowed us to analyse, for the first time, the 
deceleration associated with drift motion opposite 
to the solar wind electric field. We showed that 
the overall deceleration experienced by SEPs is 
due to contributions from adiabatic deceleration 
(which requires scattering) and drift-induced de¬ 
celeration (which is present also in the absence of 
scattering). 

We found that energy changes associated with 
the combined processes of drift-induced and adia¬ 
batic deceleration are very large: for a scattering 
mean free path of A=1 AU, AK/Kq at 100 hours 
after injection ranges between —35% and —90% 
for particles with A' 0 =l MeV, between —30% and 
—70% for A 0 =10 MeV and between —20% and 
-55% for A 0 =100 MeV. 

It is interesting to note that while the magni¬ 
tude of drift, as measured by the change in colat¬ 
itude, is very small for protons at ~1 MeV, its ef¬ 
fect on particle energy is large, meaning that drift 
effects cannot be ignored even at these low ener¬ 
gies. 

With reference to the very large values of de¬ 
celeration obtained from our simulations for par¬ 
ticles of initial energy of 1 MeV, e.g. values close 
to 100%, the test particle approximation breaks 
down once the particle energy has reached values 
close enough to that of solar wind particles, so that 
a different modelling approach is required. 

Deceleration is fast during the first ~20 hours 
following injection and continues at a slower rate 
in the subsequent hours. The fact that it is tak¬ 


ing place also for scatter-free propagation means 
that it cannot be ascribed to standard adiabatic 
deceleration only. In addition, the simulations 
with scattering taking place in the fixed frame 
(i.e. with scattering present but adiabatic decel¬ 
eration ‘switched off’) demonstrate the large mag¬ 
nitude of drift-induced deceleration. 


The magnitude of the constant £0 appearing in 
the expression for the solar wind electric poten¬ 
tial (Eq. d5j) ) determines the particle energy range 
over which the effect of drift-induced deceleration 
is most prominent. It appears that for protons this 
is the range between ~1 and a few hundred MeV, 
which is the typical SEP energy range. At proton 
energies larger than 10 GeV the energy change 
becomes negligible, therefore its importance for 
galactic cosmic rays is for modulated cosmic rays 
in the solar wind with energies in the ~1 MeV - 
1 GeV range. 

We also emphasize that drift additional to that 
due to the large scale Parker spiral field and any 
other process which is able to produce transport 
in the latitudinal direction will cause changes in 
kinetic energy. 


Is drift-induced deceleration inclu ded within 
the standard transport equations ( Skillind 1971 


Ruffolo 1995 ) used in SEP modelling, data anal¬ 
ysis and Space Weather forecasting? Answering 
this question requires analysis of th e deri vation of 
the focussed transport equation ( Skilling! Il97lh . 
This equation describes the evolution of the dis¬ 
tribution function of the particles’ guiding cen¬ 
tres and is obtained by averaging the Vlas ov 
eq uatio n over gyrophase (e.g. IZankl ( 2014) ). 
le R oux & Webb ( 2009 ) showed that the focussed 


transport equation is equivalen t to the guiding 
centre kinetic equation ( Kulsrudlll983 ). The latter 
is derived by means of the standard assumptions 
of first order guiding centre theory, including the 
assumption that the grad B and curvature drift 
velocities are much s maller than the E x B dri ft ve¬ 
locity. As shown by le Roux fe Webb ( 2009h . the 
guiding centre kinetic equation includes a con¬ 
tribution from grad B and curvature drift in its 
energy change term, while the latter drifts are 
not included in the spatial convection term. This 
means that the energy change associated with 
drift is taken into account in an approximate way, 
with the deceleration rate being calculated along 
the ‘undrifted’ orbit, at constant colatitude. To 



































our kno wledge, all current SEP transport mod - 
els (e.g. Zhang et, al. ( 20091) . Mason et al. (2012)) 
make use of this approximation, which we expect 
to be fairly accurate for low energy SEPs but less 
accurate for high energy ones, due to their sub¬ 
stantial latitudinal drift. A more precise solution 
would be obtained by solving a transport equa¬ 
tion that includes the grad B and curvature drifts 
within the spati al co nvection term, i.e. by solving 
eqs. (30)-(34) of lZhand (l2006h . 


It should also be noted that in Paper II we 
showed that the grad B and curvature drift ve¬ 
locities in the Parker spiral can be of the same 
order of magnitude as the E x B drift velocity, and 
in some cases even larger. Therefore the standard 
guiding centre theory assumption that the ExB 
drift is the dominant one may not be valid for the 
case of heliospheric propagation. 


We conclude that the extent to which the drift- 
induced deceleration observed in our test particle 
simulations is accounted for by SEP models based 
on the focussed transport equation needs to be 
evaluated, and this will be the subject of future 
work. 


Strong SEP deceleration such as the one we 
reported, has an impact on the analysis of mea¬ 
surements, because it introduces a scenario where 
particles detected in a given energy channel at 1 
AU actually left the Sun ha ving a much higher 
energy, as was emphasized bv lMason et al. ( 2012 ) 
for heavy ions in the ~MeV/nucleon range. Our 
results show that this may be taking place for pro¬ 
tons in the 10-100 MeV range too. 


Overall, our results demonstrate that drift- 
induced deceleration is an important process in 
SEP propagation and needs accurate inclusion 
within SEP modelling and data analysis. 
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